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Abstract 

We study the power spectrum of super-Hubble fluctuations of an inflaton-like scalar field, the 
“system”, coupled to another scalar field, the “environment” during de Sitter inflation. We ob¬ 
tain the reduced density matrix for the inflaton fluctuations by integrating out the environmental 
degrees of freedom. These are considered to be massless and conformally coupled to gravity as a 
proxy to describe degrees of freedom that remain sub-Hubble all throughout inflation. The time 
evolution of the density matrix is described by a quantum master equation, which describes the 
decay of the vacuum state, the production of particles and correlated pairs and quantum entangle¬ 
ment between super and sub-Hubble degrees of freedom. The quantum master equation provides a 
non-perturbative resummation of secular terms from self-energy (loop) corrections to the inflaton 
fluctuations. In the case studied here these are Sudakov-type double logarithms which result in 
the decay of the power spectrum of inflaton fluctuations upon horizon crossing with a concomitant 
violation of scale invariance. The reduced density matrix and its quantum master equation fur¬ 
nish a powerful non-perturbative framework to study the effective field theory of long wavelength 
fluctuations by tracing short wavelength degrees of freedom. 


*Electronic address: boyan@pitt.edu 


1 




I. INTRODUCTION 


Cosmology is entering the precision era, measnrements of the cosmic microwave back- 
gronnd anisotropies by the WMAP[I| and PLANCK missions have provided conhrmation 
of a nearly scale invariant power spectrnm of adiabatic pertnrbations, one of the main pre¬ 
dictions of inflation. These measnrements snpport the main paradigm of inflation based on a 
scalar held slowly rolling down a potential landscape leading to a nearly de Sitter inhationary 
stage and the generation of pertnrbations from the qnantum hnctnations that are amplihed 
when their wavelengths cross the Hnbble radius. This simple and compelling picture invites 
deeper scrutiny at a fundamental level. Most models of inhation rely on either a single scalar 
held or several scalar helds (which typically yield isocurvature or entropy perturbations) but 
no interactions between the inhaton held and the large number of degrees of freedom of the 
standard model of particle physics (and beyond). Excitation of the degrees of freedom that 


would populate a radiation dominated era 
of inhation during a period of “reheating' 


ost inhation is assumed to happen at the end 
Sl, however, this necessarily implies a coupling 


between the inhaton and the degrees of freedom that describe the physics of the radiation 
dominated era. Such interaction between the inhaton and other (fermionic, scalar) degrees 
of freedom cannot just switch-on at the end of inhation, and on physical grounds should be 
expected to be present even during the inhationary stage. Interactions of quantum helds in 
de Sitter (or nearly de Sitter) space-time have been the focus of several studies jbl-l^ which 
show strong infrared and secular ehects and the possibility that the vacuum state in de 
Sitter space time may be unstable towards decay [2^. Non-Gaussianity is a manifestation 
of self-interactions of curvature perturbations and could leave an observable imprint on the 
cosmic microwave background, although it is argued to be small in single held slow roll 
inhat ionary mo dels (2^291 . 

Interactions with heavy helds with masses M ^ H with H the Hubble parameter during 


inhation have been treated in terms of ehective held theory descriptions |30N33l| mainly by 


neglecting kinetic terms and correlations, ehectively treating the heavy degrees of freedom 
as auxiliary helds that can be “inte grat ed out” at tree level, or including correlations of the 
heavy helds in powers of H/M -C ijSdj. 

In a non-equilibrium situation as is the case with cosmological expansion, integrating out 
short wavelength degrees of freedom leads to the ehective held theory description in terms of 
a reduced density matrix for long wavelength hnctnations. Such a description is, fundamen¬ 
tally, akin to a Wilsonian approach to an ehective held theory}^ by coarse graining short 
distance degrees of freedom. At the level of a non-equilibrium ehective action, the study 
of the ehects of tracing out degrees of freedom was pioneered with the study of quantum 


Brownian motion[36H39[, the degrees of freedom of interest are considered to be the “system” 
whereas those that are integrated out (traced over) are the “bath” or “environment”. The 
ehects of the bath or environment are manifest in the non-equilibrium ehective action via 
an influence action which is in general non-local and describes dissipative processes. This 
inhuence action is determined by the correlation functions of the environmental degrees of 
freedom, and determines the time evolution of the reduced density matrix. An alternative 
but equivalent description of the time evolution of the reduced density matrix is the quantum 


master equation\4:(t l41l| which includes the ehects of coupling to the environmental degrees 


of freedom via their quantum mechanical correlations. 

A generic quantum master equation approach for a reduced density matrix describing 
cosmological perturbations has been advocated in ref. ^ in terms of local correlations of 
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environmental degrees of freedom. 


In ref. the equivalence between the influence functional and the quantum master 


equation in Minkowski space-time was established, and shown that they provide a non- 
perturbative resummation of self-energy diagrams directly in real time providing an effective 
held theory description of non-equilibrium phenomena. 


Motivations and goals: In this article we consider an inhaton-like scalar held as “the 
system” in interaction with other helds considered as the “environment” with the goal of 
studying the inhuence of sub-Hubble degrees of freedom of the “environment” upon the 
power spectrum of super-Hubble huctuations of the “system” during de Sitter inhation. 
The quantum huctuations of the inhaton-like scalar (the system) are amplihed and become 
classical when their wavelengths become larger than the Hubble radius during inhation and 
in the non-interacting theory their power spectrum becomes nearly scale invariant. We 
study the inhuence of the environmental helds by obtaining the reduced density matrix 
for the inhaton-like scalar held, and the quantum master equation that describes its time 
evolution by consistently tracing over the environmental degrees of freedom. Our goal is to 
obtain the corrections to the power spectrum of super-Hubble huctuations of the inhaton- 
like scalar from the interaction with degrees of freedom whose quantum huctuations remain 
sub-Hubble all throughout the inhationary stage. For this purpose we consider the system 
to be a minimally coupled scalar (inhaton-like) coupled to a massless conformally coupled 
scalar held-the “environment”- that serves as a proxy for helds (including fermionic helds) 
whose quantum correlations do not become amplihed for super-Hubble wavelengths. We 
also comment on the case of minimally coupled environmental helds. The mode functions 
of the quantum environmental helds ehectively describe quantum huctuations that remain 
sub-Hubble all throughout inhation. 


Brief summary of results: We obtain the reduced density matrix of the inhaton- 
like scalar held (the system) and its quantum master equation, by tracing out sub-Hubble 
degrees of freedom (the “environment”) up to second order in the coupling. We consider 
the case in which both the “system” and the “environment” are in their Bunch-Davies 
vacuum states at the beginning of the inhationary stage. A perturbative analysis of the 
quantum master equation explicitly shows the decay of the vacuum state and the production 
of single particles as well as correlated pairs which lead to quantum entanglement between 
the inhaton huctuations and those of the environmental helds. The full solution of the 
quantum master equation provides a non-perturbative resummation of self-energy diagrams 
determined by the correlation functions (loops) of the environmental degrees of freedom. 
From the quantum master equation we obtain the equations of motion for super-Hubble 
correlations of the inhaton held from which we extract the power spectrum. Its solution 
provides a non-perturbative resummation of secular Sudakov-type double logarithms from 
the inhaton self energy and yields the corrections to the power spectrum of super-Hubble 
huctuations. These indicate the decay of the power spectrum after “horizon crossing” and 
violation of scale invariance even when the power spectrum in absence of interactions is scale 
invariant. 


II. THE MODEL: 

We consider a spatially hat Friedmann-Robertson-Walker (FRW) cosmological space-time 
and two interacting scalar helds 0, ip although the methods and broad conclusions will be 
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more general. The field 0 is an infiaton-like scalar field minimally coupled to gravity, this is 
the “system”, and the field 99 is the “environment” as discussed below it will be chosen to 
be a massless, conformally coupled scalar field and it will be traced out of the total density 
matrix to yield the reduced density matrix for 0 . 

In comoving coordinates, the action is given by 


S = J (fx dt a^{t) I ^0^ - ^ (^M| + 


with 


7^ 



( 2 . 2 ) 


being the Ricci scalar, .^ = 0,1/6 correspond to minimal and conformal coupling respectively. 
The interaction has been normal-ordered 


: 45^ := 4.^ - 


(2.3) 


where the brackets ((■ ■ ■)) refer to the expectation value in the initial density matrix (see 
below). 

In the case of de Sitter space time with a{t) = it is convenient to pass to conformal 
time f] = with 

a{t{ri)) = 


Hr] ’ 


and introduce a conformal rescaling of the fields 


X{x,v) 

= 77773T ; ‘fyy) = 


a{t{7])) 

After discarding surface terms the action becomes 


a(t(T/)) 


(2.4) 


(2.5) 


^ = / dy dr] [x'" - (vxy - Mliv) x"+y" - mr - m^v) y 


Hr] 


X ■ 


with primes denoting derivatives with respect to conformal time r] and 


= 


H2 +12(^^x,b g 




where for consistency of notation we have called 


(2.6) 

(2.7) 

(2.8) 


respectively. 
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Since A has dimensions of mass we will consider the weak coupling case with \/H 1. 
In the non-interacting case A = 0 the Heisenberg equations of motion for the spatial Fourier 
modes of wavevector k for the conformally rescaled fields are 


x'M + 


^ a\ X 47j 


+ 


7]^ 

1 

?7^ 




Xfc(h) = 0 
= 0 


where 


9 /M2 


“ 4 Vhf2 


+ 12^a) ; a = x,^- 

The Heisenberg helds are expanded in a comoving volume V as 

x{x,v) = -^'^[bq9{q,v)+ bl^9*{(l,v) 


+ alqU*{q,ri) 


W 

We choose Bunch-Davies conditions in both helds, namely 

%|0)x ~ 0 j ~ 0 


Ak-x 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 


and 

9(9, V) = (2.15) 

u{k,ri) = i y/—nriHl^J{—kri). (2-16) 

These conditions may be generalized to non-Bunch-Davies, but here we consider this simpler 
case to highlight the main physical consequences. 

The X held is considered to be minimally coupled, ■Cx ~ 0 nearly massless with 
M-^/H 1, from which it follows that as —qrj —)■ 0 


9 A 9 ) « 1/h- 


(2.17) 


This behavior in the super-Hubble limit will lead to strong secular contributions in the long 
time limit. 

The time evolution of a density matrix initially prepared at time ? 7 o is given by 

pA) = UA,Vo)p{vo)U-\ri,rio ), (2.18) 


where Tr[p(? 7 o)] = 1 and U{r],r]Q) is the unitary time evolution of the full theory, it obeys 


i-^U{r],7]o) = H{7])U{7],7 ]o) ; H(r/o,r/o) = l 


(2.19) 


where H{r]) is the total Hamiltonian. Writing the total Hamiltonian in terms of the free and 
interaction Hamiltonians as H{r]) = HqA) + Hi{r]) it is convenient to pass to the interaction 
picture introducing the unitary time evolution operator of the free theory UoA,r]o) obeying 


i-^UoA, r]o) = Ho{ri) Uo{t], rjo) ; Uo{t]o, t]o) = 1 

drj 


( 2 . 20 ) 
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and define the density matrix in the interaction pictnre 

Pi{v) = Uo\v,Vo)p{v)Uo{r],r]o) = Ui{r],r]o)p{r]o)Uf\r],r]o) ( 2 . 21 ) 

where Ui{'q,r]Q) is the nnitary time evolntion operator in the interaction pictnre obeying 

i^Ui{r],r]o) = Hj{r])Ui{r],r]o) ; f//(?7o, ho) = 1 (2.22) 

where 

Hi{p) = Uq^t], r]o)Hi{r])Uo{r], r/o) • (2.23) 

From fl2.ip we hnd 

Hiiv) = J d^xxix,v) : V’^(x,h) : (2.24) 

and the y, '0 fields are in the free field Heisenberg representation fl2.12ll2.13p . Normal ordering 
is defined as 

: 7]) := 'ijj'^{x, h) - Tt[iP^{x, h)p/(ho)] (2.25) 

III. QUANTUM MASTER EQUATION 

The steps leading to the qnantnm master eg nation np to second order in the conpling are 
given in detail in Minkowski space time in ref.j^. In this reference the eqnivalence between 
integrating ont the heavy (or short wavelength) degrees of freedom in the path integral 
representation and the qnantnm master eqnation is established np to second order in the 
conpling. This eqnivalence and the relation to a stochastic description translate directly 
to the case of an FRW cosmology, these more formal aspects are relegated to a companion 
article ji^. In this article we obtain the qnantnm master eqnation directly and apply it to 
nnderstand several physical conseqnences. 

The time evolntion of the density matrix in the interaction pictnre is given by 


p'liv) = -'i[Hi{v),pi{v)] 

whose formal sointion is 

(3.1) 

rv 

Piiv) = PiiVo) - i / dp', 

this sointion is inserted back into fj3.ip leading to the iterative eqnation 

(3.2) 

p'liv) = -i[Hiiv), PiiVo)] - ^ 

f [Hiiv),[Hiiv'),Piiv')]]dv' ■ 

Vo 

(3.3) 

The next steps leading to the qnantnm master eqnation rely on various 
discussed in detail in 43 . The first is factorization, namely 

approximations. 

Piiv) = Pixiv) ® PipiVo) 

(3.4) 


this is an nbiqnitons approximation jdO. 1^ . In ref.ji^ it is shown that this approximation 
resnits from obtaining the non-eqnilibrinm effective action in the path integral Schwinger- 
Keldysh formnlation in a consistent cnmnlant expansion of the trace over the environmental 
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degrees of freedom when the initial density matrix is factorized. These aspects translate 
directly to an FRW cosmology in conformal time as discussed in a companion article ^ . In 
this article we assume (as is common in the literature) that the initial density matrix at the 
beginning of inflation is of the factorized form, it is clearly interesting to consider the case 
of initial correlations, which is postponed to further study. 

Taking the trace of pi{r]) over the V’ degrees of freedom yields the reduced density matrix 


PrM = Tq,,) p,(r]) . 


(3.5) 


The normal ordering fl2.25p entails that the hrst term in fl3.3p vanishes upon taking the trace 
over 'll) leading to 

y2 rv dr)' f f \ 

pr('d) = ~rip j j ~ 

+ pr{v')x{y, v')x{x, v) G<{x - y, r], r]') - x(x, v)pr{v')x{y, v') G<{x - y, r], r]') 

- x{y,v')Pr{v')x{x,v)G>{x-y,v,v')'>, (3.6) 


where 


G'-'li- ff,= Tr[: >?') : Prtivo)] , 

G^(x-y, 'ii 'if) =Tr|: ■■ : /3/*(i)o)|. 


(3.7) 

(3.8) 


In writing the correlation functions as functions of T — y we used spatial translational 
invariance in a spatially flat FRW space-time which allows us to write. 


1 


{x-y, V,v') = y^lC>[p,r], r]' 


^-ip-{x-y) 


1 


; G< {x-y, ri, V') = yYl 


„-ip-{x-y) 


(3.9) 


where V is the quantization volume. 

A this stage, the second Markov approximation is invoked: taking pr{p') —t Pr{v)- This 


is justihed in weak coupling: since dpr{p)/dp cx 1, an integration by parts |43 

and neglecting contributions of 0{\'^) on the right hand side of fl3.6p leads to the Markov 
approximation. The main arguments showing the validity of this approximation for weak 
coupling are available in ref. 43] . 

In order to obtain the correlation functions we need to specify pip{po), we consider 
the simple case 

PitpiVo) = |0)i/,p(0| (3.10) 

where |0)^ is the Bunch-Davies vacuum for the "0 helds and —pQ is taken to be the beginning 
of the (nearly) de Sitter inflationary stage, although a generalization is straightforward. In 
this case we hnd 


IC>[q]p,p'] = K[q]p,p'] = 2 j 
/C< (g; p, p') = K> [q- p', p] = K* [g; p, p'] , 


u{k,p)u*{k,p')u{p,p)u*{p,p') ; p = \k + (^ 

(3.11) 


where u{k,p) etc, are the mode functions fl2.16p . The Feynman diagrams for the correlation 
functions G^,G^ are shown in £g.([T]). 
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FIG. 1: The correlation functions (x — y, rj, rj'), {x — y, y, y'). 


Using the expansion fl2.12|) we find the general form of the quantum master equation up 
to second order in the interaction fl2.24l) 


p'ri.y) = 


A2 p dy' 


H^vJvo y 


— 1 Xqiv)x-qiv')pr{v) K{q-,y,y') + priv)X-qiv')Xq{v) K*{q]y,y') 


Xq{v)pr{v)X-q{v') K* [q] y, y') - X-q{v')pr{v)Xq{v) K{q-, y, y') 


(3.12) 


where 

Xq{y) = b^giq, v) + v) ; x-giv) = xliv) ■ (3.13) 

To obtain expectation values of operators in general one has to solve the quantum master 
equation and hnd Pr{y)- However, it is most useful to obtain the evolution equations for 
operators that do not evolve in time in the interaction picture. Consider one such operator 
O which is constant in the interaction picture, then with {0){y) = TiOpriy) follows that 

(3.14) 

It proves illuminating to write the quantum master equation in terms of the operators 
b]j,bq, this yields simpler expressions for the matrix elements of the reduced density matrix 
in the Fock basis of quanta associated with these operators, in the case under consideration 
these are Bunch-Davies Fock states. We hnd 


Priv) = 


dy' 


H^y 


'VO 


7 1 y)K[q, y, y'] bjjb^priy) - b^priy) bl 


9{(i,v)9*{(i,y,y'] priv - b^pr{y) bi 
giq, y)9*{(l, y')K[q] y, y'] b^lpriy) - blpriy) b^ 
9{(l,y')9*{(l,y)K*[q]y,y'] Pr{y) bpl-bipr{y) b^ 
gig, y)9ig, y')K[q] y, y'] b^-^priy) - b-^priy) 6g- 
9ig,y)9ig,y')K*[q]y,y'] priv)bp-^-b_^priy)b^ 
9*ig, y)9*ig, y')K[q-, y, y'] by)Kpriy) - bKpriy) 6t 

9*ig, y)9*ig, y')K*[q; y, y'] \priy) b^qblff - blffPriy) bl 


(3.15) 
























Introducing the expectation values 


Nqiv) = {b^bq) = '^^(b^bqPriv)) i = (b^-q) = (b^-q Pr{v)) (3-16) 


we find 


KM = + +-Sf’W 

KM = ia,MKM+2il3’(.r))N,(r)) + S!K'’M, 


where 


Iqiv) = 

si^\v) = 
I3qiv) = 
Sf\p) = 


aq{ri) = 


2 d7]' 


4A 

Ir^ 

A2 
H^t] 

2A2 n d7]' 


'vo 

n di r 


jKi[q]Tq,Tq']\m g{q,v')g*{q,v) 


(3.17) 

(3.18) 

(3.19) 


Mo 


g' 


g{g^ v)g*{g, v')K[q-, g, g'] + g{q, g')g*{q, g)K*[q; g, g'] (3.20) 


H^g 


m 

2A2 n di 


yKi[q]g,g'] ^g{q,g')g{q,g) 
YK[q-,g,gy*{q,g')g*{q,g) 


H^vJqo V 
^ V Jqo V 


K[q;g,g'] = KR[q;g,g'] + iKj[q;g,g'[ 


(3.21) 

(3.22) 

(3.23) 

(3.24) 


Once the kernel K[q]g,g'] is found the above equations can be integrated. The in¬ 
homogeneous source terms Sq^\g)-, Sq^\g) are noteworthy, if these vanish, the vacuum 
Nq = 0, Mg = 0 would remain a fixed point of the dynamics, therefore these source terms 
indicate particle production and the production of correlated pairs of particles. Since these 
terms are independent of the initial number of particles and emerge even when the initial 
density matrix corresponds to the vacuum pure state, they can be understood perturbatively. 


IV. PERTURBATIVE INTERPRETATION: PARTICLE PRODUCTION, AND 
ENTANGLEMENT 

Before proceeding to analyze the full quantum master equation, it proves illuminating to 
understand the different contributions from the point of view of a perturbative expansion in 
the coupling. This analysis yields important insight into the correlations between the system 
and the enviroment and how the quantum master equation provides a non-perturbative 
resummation of self-energy contributions. 

Consider that the initial density matrix describes the pure vacuum state pi{go) = |0)(0| 
where |0) = |0)^ |0),/,, up to second order in the interaction we find the state 

|1,(,,)> = |0> + |>1.C)(,,)) + |>1.P>(,,)> + ... (4,1) 

where 

I*'"!’))) = (-1) rH,M)dq’\0), (4.2) 

p rvi rv 

\'^b^\g)) = (-if Hj(gi) Hi(g2)dgidg2\0) = (-i) Hi(gi)\'^^^\gi)) dgi(A.3) 

-^rio Jr}o 
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Up to second order the reduced density matrix up is given by 


Pr(h) = Tr^l^(h))(^(^)l = Tr^ 


|o)(o| + 


• 4 ) 

There is no cross term | 0 )(\I/^^^(r 7 )| since the state |\I/(^)(r 7 )) has vanishing overlap with |0)^ 
because of the normal ordering in the interaction Hamiltonian (see below). 

We hnd 


f —9*{(l,Vi)u*{k,r]i)u*{p,r]i)\l^)^\lj:,lj^)^ ; p=-k-q, 
HVV . . Jr,o Vi 

(4.5) 

where |l^)x are single y particle states and similarly for the kets with -0 particles. The state 
fl4.5p is an entangled quantum state of the system (y) and environmental (-0) helds. This 
state describes particle production from the vacuum state by the interaction, depicted in 
£g.(l2]). In Minkowski space time these are virtual processes as they do not conserve energy, 
however in an expanding cosmology these processes are available as “real” because there is 
no time-like Killing vectorj^, [^. In particular the mode functions g{q,rii) feature a growing 
component when q crosses the Hubble radius g{q,p) — leading to secular (growing) 

contributions to the time integrals as 77 —)■ 0 for z/^ 3/2 for a minimally coupled scalar 

held. This feature will become important below when we discuss the power spectrum. 


q 

X 







FIG. 2: The state |'I'd)(ry)): production of correlated x ^md V’ particles from the vacuum state. 


In second order the state |T*^^^(? 7 )) features several contributions obtained by applying 
the interaction Hamiltonian to |Td)( 7 y)) as per the second equality in fl4.3l) . However, only 
two of these contribute to Pr{g) to second order, these are: I) annihilate all particles in 
| 4 /d)( 7 ^)) returning to the full vacuum state |0), H) create another y particle and annihilate 
(both) 'll) particles returning to the -0 vacuum state but to a two particle state of the y held. 
Only these two contributions feature an overlap with the vacuum state necessary for the 
third and fourth terms in fl4.4p . These two states are given respectively by 

= [ — [ — 5^£/*(?,hi)^(?,h2)^[g;h2,hi]|0)x|0)p (4.6) 

-n JrjQ P2 Jrio Vl - 

=-^ [ —[ —J23*(9,Vi)9*{(l,V2)K[q-,p2,Vi]mA-q)x\^)p- ( 4 - 7 ) 

Jrjo V2 Jr^o 9l - 
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Ip tp 

(a) (b) 


FIG. 3: The states in that contribute to Pr{r])'- (a) = ; (b) = . 


and their associated Feynman diagrams are shown in hg.([3]) 

The state n describes the production of correlated particle pairs of the held y 

out of the initial vacuum state. Inserting these results into fl4.4p and carrying out the trace 
over the ifb degrees of freedom, we hnd up to second order 


Priv) = | 0 )( 0 | + ^ 


(g; p) 11,-) (1,H + Gf (g, p) \ 1.^) (01 + (Cf ^ (g, p))*] 0) 1 


(2), 


9 ’ -^9 


+ C'y(g,.,)|0)(0| + (Cf'(,,.,))10){0| 


(4.8) 


where now all the kets correspond to y particle states (we suppressed the label y to simplify 
notation). The coefficients are given by 


C^^\q;v) = ^ [ — [ —9iq,V2)g*{q,Vi)K[q;p2,Pi]>0 (4.9) 

Jr,o V2 Jno Vi 

Cf\q,p) = f ^ f —g*{q,pi)g*{q,p 2 )K[q;p 2 ,pi] (4.10) 

^2 Jr,o Vl 

Co\q^v) = -4t [ — [ —9{q,V2)g*{q,Pi)K[q;p2,pi] (4.11) 

J r}o V 2 J r]o Vl 


The coefficient C^^\q]p) can be written in a more illuminating manner by introducing 
0 (h 2 ~ hi) + 0(hi ~ V 2 ) = 1 in the integrals, in the term with Q{pi — P 2 ) relabel pi O p 2 
and use the property K[q-,pi,p 2 ] = iF[g; 172 , hi]* fo hnd the relation 

C'(^)(g;h) = -(cl^\q,p) + iC^^\q,p))*^ . (4.12) 

This identity conhrms unitarity in the total time evolution, since 


TrxPr(h) = ^ + S [C^^\q',v)+ Co\q,p) + {C^^\q,p)T 

9 


1 = Tip(r;o). 


(4.13) 


Several aspects stem from the above results: 

• i): with blb^= Ng it follows that 

TiNgPrip) = {Ng)ip) = C^^\q;p) , (4.14) 
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the relation fl4.12p leads to 


KM = 




H'^ r] 


'Vo 


dr]i 
Vi L' 


d^k,V,Vi]+9*i(l^v)9i(l,Vi) d<*k,V^Vi] , (4-15) 


this is precisely the source term Sq^\ri) fl3.20p in the rate equation fl3.17p . and makes 
manifest the production of x particles from the vacuum state, namely the decay of the 
vacuum during the time evolution. Whereas in Minkowski space time these are virtual 
processes leading to the wave function renormalization of the vacuum, the lack of a 
time-like Killing vector as a consequence of the cosmological expansion makes these 
processes to contribute in the long time limit 10 . 


The super-Hubble limit of the mode functions (in the nearly massless case) 9 ( 9 , 9 ) oc 
1/9 fl2.17p leads to secular growth of particle production and vacuum decay in per¬ 
turbation theory. These secular contributions are effectively re-summed by the full 
quantum master equation. This statement will become clear below. 


ii): The probability of the x vacuum state in the reduced density matrix is 

Z{ri) = (0b,('i,)|0) = 1 -^K(ri ), (4.16) 

which is interpreted as the vacuum wave function renormalization, namely the prob¬ 
ability of hnding the “bare” vacuum state in the full state time evolved from the 
vacuum. This result is in agreement with that found in ref. 21j within the Wigner- 
Weisskopf approximation, if the number of particles does not saturate in time, the 
Bunch-Davies vacuum is unstable towards particle production consistent with the 
results of refs. UM Hi 


The expectation values M*{vi) = describe the pro¬ 


duction of correlated pairs. 


• iii): The contribution |Td)(r^)) given by fl4.5p and in the trace fl4.4l) clearly shows that 
the degrees of freedom of the system x are entangled with those of the environment 'ip. 
In particular for q -C —1/?7 and k ^ —1/?7, there is quantum entanglement between 
super-Hubble modes of the system and sub-Hubble modes of the environment. This is 
in agreement with the super-sub-horizon entanglement discussed in ref. [l^ . 


• iv): Although the last two terms in fl4.8p can be combined, they arise from 
|4/*^^^(?7))/(0| and | 0 )/(T*^^^(? 7 )| respectively, then writing the various contributions in 
04.81) in terms of bt, b^ we hnd 


Pr{9) = Pr {90) + 


1 L 


C^^\q; 9 ) bipri9o) b^+C^^\q, 9 ) + 


(Cf\q, 9))* pr{9o) feq + C^^\q, P) bqb\pr(9o) + (C^o\q, 9))* PriPo) 4^9 


(4.17) 


Upon using the relation 04.12|) and taking the derivative with respect to 9 one recog¬ 
nizes by inspection the similar terms in 03.15P when taking Pr(h) PriPo) = |0)(0| on 
the right hand side of 03.15P to lowest order in perturbation theory. 


Therefore it becomes evident that 03.15P along with 03.17113.18p furnish a non-perturbative 
resummation for the time evolution, similar to quantum kinetic equations as discussed in 
ref. 43|] in Minkowski space-time. 


12 








































V. TRACING OUT SUB-HUBBLE MODES: 


The discussion in the previous section has been quite general and does not specify in 
detail either the mass or the coupling to gravity of the environmental degrees of freedom -0, 
namely the value of in fl2.7p . 

Our main goal in this article is to understand the effective evolution of super-Hubble 
fluctuations of the system held y (inhaton) upon integrating out (tracing over) huctuations 
of the environmental held with wavelengths that remain sub-Hubble all throughout inhation. 
In other words, we seek to obtain an ehective held theory for long wavelength modes, those 
that would be of cosmological relevance today, tracing out short wavelength modes of the 
environmental helds. For physical wavelengths that are much smaller than the Hubble radius 
—kr] 3> 1, the Bunch-Davis mode functions of the environmental helds 

_ y Q-ikri 

u{k, rj) —kr] 3> 1 

If the environmental helds are scalar helds minimally coupled to gravity, their quantum 
huctuations feature a growing mode that is amplihed upon becoming super-Hubble. If the 
X is the source of (adiabatic) perturbations a minimally coupled environmental scalar held 
ip would yield isocurvature (or entropy) perturbations, which are constrained by CMB ob¬ 
servations. This motivates us to consider the bosonic held ip to be massless and conformally 
coupled to gravity, namely = 1/6 in which case = 1/2 and 



u{k,ri) = 


Q-ikr] 


(6.2) 


for all k,r]. This is also the case if the bosonic held y couples to a fermionic held with 
mass mj -C i/, the mode functions for the fermionic helds are those of Minkowski space 
time but in terms of conformal time[l^, 45, ^|. Therefore, choosing -0 to be a massless 
conformally coupled scalar held is a proxy for integrating out (tracing over) sub-Hubble 
degrees of freedom, leaving an ehective action for the degrees of freedom that become super- 
Hubble during inhation. With the mode functions fl5.2l) it is now straightforward to obtain 


K[q,r],r]'] 


i Q-igiv-v') 
Svr^ {rj — rj' — is) 


£ ^ 0 + 


(5.3) 


where e is a short-distance cutoh that regulates the momentum integral in fid.lip . This result 
can also be conhrmed from the operator product expansion of the composite operator : ip"^ : 
in Minkowski space-time since the mode functions are the hat space time plane waves and 
the normal ordering is in the vacuum state. If, instead, we considered the environmental 
held [ip) to be minimally coupled, there is an additional contribution to K\m V. v'] from 
modes that become super-Hubble, this is discussed in the second reference in|22|] and will 
not be considered further here. 

Therefore K[q,r],r]'] features a local and a non-local contribution. 


K[q,T],T]'] 


871 ^ 


V 


- 1 - 
-T] — 7]'. 




(5.4) 


where V stands for the principal part. The result fl5.4p is important, the correlation functions 
of environmental helds are often taken to be local in time, namely oc 6[r] — p') as in the 
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second term in fl5.4p . however even for fields with sub-Hubble fluctuations there is a non¬ 
local contribution which describes a long range memory. As it will be shown below, this 
non-local contribution is of paramount importance in the quantum master equation, it is also 
present if the environmental helds are minimally coupled to gravity (see second reference 


m 


221 ). 


It is convenient to write the non-local term in fl5.4p as 


V 


r 1 1 

p — p' ^ 1 

{p — p'Y -I- 

-P — p'. 

{p — p'Y + 2 dp' 

{-pY 


(6.6) 


In fl5.5l) we have introduced an arbitrary scale {—rf) to render the argument of the logarithm 
dimensionless and acts as a subtraction or renormalization scale just as in the usual renor¬ 
malization program in Minkowski space-time. It will be judiciously chosen below. We can 
now input fl5.4p with fIS.Sp in fl3.12p and separate the local contribution from the S{T]—r]') and 
the non-local contribution from the principal part to the quantum master equation fl3.12p . 
The non-local contributions feature the integrals 


[ K[q;T],T]']x-q{v')^ ; [ K*[q;T],T]']x-q{v') ^ , 

'VO ^ J Vo ^ 


(6.6) 


the form fIS.Sp allows to extract the divergence from the non-local contribution upon inte¬ 
gration by parts. 

At this point we choose the arbitrary renormalization scale ^ = tjq with —rjo very large 
(the beginning of the inflationary stage) so that the wave vectors q of interest are sub-Hubble 
at this time, namely —qr]o ^ 1. When evaluated at rjo the lower limit of the integrals (I5.6p , 
and for —rj —)■ 0 the logarithm in (15.51) yields a contribution ~ T^/r^o which, when combined 
with the large —qr] behavior of g{q,rj) yield a contribution of order /[H qrjo)'^ 1 to the 
quantum master equation. Therefore, choosing r] = tjq in fl5.5p allows to neglect the contri¬ 
bution from lower limit in the integration by parts. The manifestation of a renormalization 
group invariance of this choice will be discussed below. 

We obtain, p'^{r]) = p'^^{r]) + pUiiv) where 


PrLiv) = 






In 


-Vo 


[x<]iv)X-q{v),Pr{v) 


E Xq{v)X-q{v)pr{v) + Pr{v)Xq{v)X-q{v) “ "^X-qiv)Pr{v)Xq{v) > , ( 5 . 7 ) 


for the local contribution, the non-local contribution is found to be 


p'rNLi.V) = i 






1 Xq{v)X-qiv)priv) “ Pr{v)X-qiv)Xq{v) 


+ Xq{v)pr{v)X-q{v) - ^-qiv)Pr{v)Xq{v) 


(5.8) 
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where 


X-M 


X-M 


'1 d 
dr]' 


'VO 


^-iq(v-v') X-qiv') 


7]' 


In 


7] — 7j 


d 

drj’ 


'VO 


^iq{v—v') ) 


In 


7] — 7] 

-Vo 


Vo J 
/ 


dr]' 

dr]'. 


(5.9) 

(5.10) 


The hrst term on the right hand side of the local contribution fl5.7p (the commutator) is 
identihed as an ultraviolet divergent mass renormalization, indeed this term can be written 
as 




In 



[xqiv)X-g{v),Pr{v) 

q 




SH{v),PrM 


(5.11) 


with 


6H{7]) 


5M\7]o) 

2 H^t]^ 


Xq{v)X-q{v) 

q 


SAiHm) 




(5.12) 


dM‘^{7]o) is an ultraviolet divergent mass renormalization of the same form as in Minkowski 
space time. This contribution is absorbed into a mass renormalization, the renormalized 
mass is 

MUvo) = M^ + 6M\r]o), (5.13) 

where we emphasized that this renormalized mass depends on the renormalization scale r]o. 
In what follows we absorb this mass renormalization in the original Lagrangian and include 
a counterterm to precisely cancel fl5.1ip thereby neglecting the first term in fl5.7p . 

We can now insert the kernel fl5.3p into the coefficient functions (13.19113.2^ in the equa¬ 
tions for Nq{r]) and Mq{r]) (13.171 IXTKp and solve them to obtain their time evolution. How¬ 
ever, while understanding their dynamical evolution merits such study on its own, neither 
Nq{r]) nor Mq{r]) are directly observable. Instead an important observable is the power 
spectrum 

v) = 4>-q{v)) = ^ ^ 2 ^ (Xg(h) X-q{v)) ■ (5.14) 

In terms of the field expansion and the expectation values (I3.16P it is given by 


V{q]v) 


rf 

27r^ 


(1 + 2N,{if}) |g(?.i()P + M,(i;)/(g,i;) + M*(i;) (9*(g,i())^ 


(5.15) 


Furthermore, only the power spectrum for wavevectors of cosmological relevance are of 
observational interest, these correspond to wavelengths that became larger than the Hubble 
radius ~ 10 e-folds before the end of inflation. The quantum fluctuations corresponding 
to these mode functions are amplified upon becoming super-Hubble and become classical, 
the mode functions feature a growing and a decaying mode and only the growing mode is 
relevant when the fluctuations become super-Hubble. This means that only a particular 
combination of Nq{j]) and Mq{7]) determined by the growing modes is relevant to the power 
spectrum. Rather than solving the coupled set of equations (13.17113.ISp it is more convenient 
to directly study the dynamical evolution of the relevant combinations. Therefore we re-cast 
the quantum master equation in terms of the degrees of freedom that have direct relevance 
to the power spectrum of super-Hubble fluctuations. 
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VI. CLASSICALIZATION AND CORRECTION TO THE POWER SPECTRUM: 


Rather than studying the dynamics of particle production and two-particle correlations 
and solving the coupled equations (13.17113.18P with the time dependent coefficients deter¬ 
mined by the kernel K[q, rj, rj'] above, we focus on understanding the impact of tracing of 
sub-Hubble modes upon physical observables, in particular the power spectrum. 

The modes of a bosonic held minimally coupled to gravity become classical when their 
physical wavelength becomes larger than the Hubble radius during an inhationary stage 
(superhorizon), this is a consequence of the fact that there is a growing mode that becomes 
amplihed and a decaying mode whose amplitude diminishes!^. This is manifest in the 
mode functions g{q,ri) fl2.15p since 

= M-qv) + iyA-(iv) ( 6 . 1 ) 


and 


Ju{z) z ^ h 



; YJz) z ^ 1) “ ( ~ ) 

\ 2 ) r(n+l) ' 



m 


TT 


( 6 . 2 ) 


In order to exhibit the classicalization of the quantum huctuations in a more direct manner, 
it proves convenient to use the real mode functions Jj,; corresponding to the decaying 
and growing modes respectively to expand the held, this is achieved by introducing the 
combinations 



These P's and Q's are canonical variables obeying the canonical commutation relations 

[qt Qd = ; [q. qi = [q?. QeI =»■ (6.6) 

Introducing the real growing and decaying mode functions 


9+ipv) = xl-^Yu^i-Qv) ; g-iQ'.v) = Ju 

we can now write the held expansion fl3.13p as 

Xqiv) = Qqg+iq-, v) + Pqg-W, v) • 


-Qv), 


( 6 . 6 ) 


(6.7) 


As discussed in refs. j47l| the classicalization of huctuations in the super-Hubble limit is 
gleaned from the relation between Xq ^-nd its conjugate momentum 

^qiv) = x'qiv) = Qq9+{q-,v) + Pqg'-iQ',v) ; [n_fc(^?), xM] = -^h,q- (6-8) 

For —qr] —)■ 0 the growing solution g+{q'., g) dominates and 


P-qiv) - Xqiv) 


9+ig'^v)\ Xqiv) 

g+ipv) 


(6.9) 
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namely the commutator between 11, y becomes much smaller than the amplitude of the 
canonical variables, a necessary condition for the classicalization of the helds. This argument 
is independent of the interaction with other helds in the theory and relies solely on the fact 
that the solutions of the Heisenberg equations of motion for a minimally coupled scalar field 
feature a growing and a decaying mode in the long time limit after the particular physical 
wavelength has become super-Hubble. This feature does not apply to either fermionic helds 
(which are never classical) nor to massless conformally coupled scalar helds (at least in 
absence of interactions). 

The power spectrum of the original (un-scaled) held is 

3 

^ (6-10) 

in the limit —qr] 1, where <f>^ is the spatial Fourier transform of the original held. From 
the scaling relation fl2.5p and fl2.4p and with Xq expanded as in fl6.7|) we hnd for —qr] —)■ 0 


V{q;r]) = 


H2 22^xr2(z/^ 


i-Qv) 


S—2i'x 


{QlQ}) 


TT 




TT 


22'^x Z/^F2(z/^ 


(Pl-Pi) (- in) 


4^'v 


( 6 . 11 ) 


With (f) being a minimally coupled scalar held with Mr -C H, where now Mr is the renor¬ 
malized mass (I5.13P (here we suppressed the renormalization scale rjo in the dehnition of the 
renormalized mass) it follows that 


Z/v — - 
^ 2 


Ml 

3H^ ’ 


( 6 . 12 ) 


therefore to leading order in Mr/H we hnd 


i/2 






(-TO)“+i (qW (-n)' 


(6.13) 


A few e-folds after the corresponding mode crosses the Hubble radius, namely —qr] -C 1 the 
second and third terms in the bracket can be safely neglected and 

_V pp2 2M^ 

V{q)r]) -qr] 0 ^( - ^ (6-14) 

When the state is the Bunch-Davies vacuum {Q^^Qq) = 1/2 and for Mr —)■ 0 one hnds the 
usual scale invariant power spectrum for a massless scalar held, 

V{q;r]) -qr] . (6.15) 

(27r)2 

In order to obtain the quantum corrections from tracing out the sub-Hubble degrees of 
freedom, the average ((■ ■ ■)) must be carried out with the reduced density matrix Pr(h)- 
The free-held Heisenberg operator fl6.7p is precisely the operator in the interaction picture 
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wherein Qcf, do not depend on time, then the average in the bracket in fl6.13l) is now 
understood with Priji) and in the limit {—qrj) —t 0 the terms with {—qr])^, {—qp)^ can be 
safely neglected, leaving only the contribution Namely we need to obtain 

2M^ , . 

'P{q-,v) = T^^\Q]PqPr{v)j ■ ( 6 . 16 ) 


Therefore it is more convenient to write Pr{p) in terms of the canonical phase space variables 
Q^, instead of the creation and annihilation operators. Furthermore, we need to extract 
the coefficient functions of Q^, P^ in the expressions X_^, in (15.9115.lOp . Since the power 
spectrum is obtained in the super-Hubble limit, namely —qp —?• 0 we can simplify the 
integrals in (15.9115. lOp by i) taking gp —)■ 0, ii) cutting off- the integrals at a time scale 
p* ~ — 1/g in the lower limits of the integrals, since the integrand is dominated by the 
late time p' ~ r; —)■ 0 region. These approximations capture the leading behavior in the 
super-Hubble limit and simplify the integrals in (15.9115.lOp i leading to 

X_^{p) = Q-qg+iq;p) + P-^g_{q;p) , (6.17) 


where 


9±{q;p) = 



9±iq] p') 

In 

-p _ fj'- 


pi 


1 -po \ 


dp' ; - g?7* ~ 1 


After renormalization of the mass by a mass counterterm to cancel the hrst term in 
(15.7p we hnd the local part of the quantum master equation to be 


(6.18) 

Prdv) 


p'rLiv) = 


+ 

+ 




^ ^ A Q ^—q Priji) P Pr{ji)Q ^—q ‘^Q—qPrip)Qq {^P+iq^P)') 


IQ'irPt'^p'^ 

q 

PqP-qPr (p) + Pr ip)PqP-q “ 2P_g'P^ ip)P^ (g_ (g; p)f 

QqP-qPrip) + prip)QqP-q “ 2P_^prip)Qq (g+(g; p)g-iq-, p)) 


P^.^Prip) + Prip)P^-<f-2Q_^prip)P^ (g+(g; g)g-(g; p)) \, (6.19) 


and the non-local contribution is 


P'rNLiv) — 


{ 9+iP v)9+iq-, p) QqQ-qPrip) “ prip)QqQ-q 


%'k‘^H‘^P 

q 

+ 9-iq]P)9-iq-,P) PqP-q prip) - Prip)PqP-q 
+ 9+iq]P)9-iq]P) QqP-qPrip) - P-qPrip)Qq 
- 9+iPV)9-iq-,P) prip) P-qQq- QqPrip)P-q 
+ 9+iq-,P)9-iq]P) PcjQ-qPrip) - Q-qPrip)Pq 


9+iq;P)9-iq;P) Prip)Q-qPq- PqPrip)Q-q 


( 6 . 20 ) 
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As mentioned above, in order to obtain = TiQ^^^pr^rj) we would need to solve the 

quantum master equation. However, since in the interaction picture Qg-is time independent, 
we obtain the equation of motion for instead, namely 

■^{QqQ-g) = , ( 6 . 21 ) 

from which we obtain by integration. 

The terms bilinear in P are proportional to {g-{q', p)Y oc in p^L and p)/p (x p 
in PrNL both are subleading in the long time limit p ^ 0 and can be neglected in the 
evaluation of We hnd that yields no contribution to fl6.2ip in the long time 

limit, hnally we obtain the remarkable result 




-r(g;p) {QqQ-q) , 


( 6 . 22 ) 


where 

>'^'9+{pp)9-{pp) 

2'k‘^H‘^P 

is determined by the non-local contribution to the quantum master equation. The above 
equation is valid at long time so that the particular wavevector q has crossed the Hubble 
radius since we have neglected the contribution from the decaying mode. Therefore the 
initial condition for the integration of fl6.22p must be set at the scale p^ ~ hence 




{QqQ-q){p) = e (Qgi5_g-)(p*). (6.24) 


We can extract the leading behavior from the integral in (I6.18p in the limit —qp —)■ 0, 
and for M "^/-C 1 setting = 3/2, namely 


9+{q\p) = ^^ ; 9-{q]p) = \ q^'‘^p^ . 
q-^/^p 3 

from which we hnd to leading and next to leading order as p —)■ 0 


(6.25) 


9+{q]p) = 


g3/2 


In 


- 1 


(6.26) 


These results yield to leading and next to leading order as p ^ 0 and M^/<C 1 


T{q-p) 


Qtt'^H'^P 


In 




(6.27) 


Carrying out the integral in (I6.24p and inserting the result into the power spectrum (I6.16p 
we hnd 


V{q]p) 


El 

27r^ 


g«(<?) H-m] 




(6.28) 


where 


a{q) 


2M|M : 

3if2 


A2 


ln[-g?7o] + 1 


(6.29) 


The result fl6.28p is noteworthy: the quantum master equation yields a non-perturbative 
resummation of secular Sudakov-type double logarithms in the long time limit. 


19 




























With M'^^rjo) given by fl5.13p and by fl5.12p it is clear that a{q) is independent 

of the renormalization scale r^o, this is a reassnring conhrmation of the consistency of the 
renormalization procedure: the effective renormalized mass changes with the renormalization 
scale, but the physical power spectrum is independent of this scale, this is a manifestation of a 
renormalization group invariance. While the q-independent term in the bracket in fl6.29p can 
be absorbed as a hnite renormalization of the mass the q-dependence of a{q) is a consequence 
of the logarithmic renormalization. It is clear that even when the non-interacting theory is 
massless and the power spectrum is scale invariant, the coupling to the sub horizon degrees of 
freedom will induce a mass via renormalization effects, the logarithmic ultraviolet divergence 
implies that the renormalized mass depends on an arbitrary renormalization scale and in 
turn this results on a g dependence of a{q). 

It remains to estimate to achieve this we would have to keep the full ex¬ 

pressions for the mode functions g±{q]Ti') inside the integrals in fl5.9ll5.10p along with the 
oscillatory factors . However for —qr]' > 1 in the integration the mode functions behave 
as cos[—qri']/y/q-,sm[—qri']/y/q and the integrands are of 0 {1) in the region of integration 
—r]o —g' ^ —V* = 1/'?- The mode functions amplify sharply for —qp' < 1 and the con¬ 
tribution from this time region dominates the integrals. Therefore the region of integration 
with —rj' > —r]^ yields perturbatively small contributions of order -C 1, consequently 

{Q^-^){ri*) = {Q^-q){Vo) + O^X^/H^) ~ 1/2 -|- 0 {X^/H‘^) where we have used that at rjQ 
the initial density matrix describes the Bunch-Davies vacuum. Therefore the hnal result for 
the power spectrum for wavevectors that cross the Hubble radius is 


V{q]r]) ~ 




(2vr) 


^a(q)\’a\-qrj\ _ 


(6.30) 


Including the second term in the bracket in fl6.29p into a hnite renormalization of Mr and 
setting the renormalized mass M^{r]o) = 0 so that in absence of coupling the power spectrum 
is scale invariant, namely 

«(?) = H-QVo] , (6.31) 

we obtain the hnal form of the power spectrum 


V{q;r]) ~ 




>2 


ln[-(j)?o]ln[-ij77]-i \n^[-qri\ 


(6.32) 


It follows from this result that the power spectrum decays when the wavevector q becomes 
“superhorizon” —qp < 1. If the wavevector q is of cosmological relevance today, it crossed 
the Hubble radius ~ 10 e-folds before r//, the end of inhation, and if inhation lasts ~ 60 
e-folds, it follows that 

ln[—gr^o] ~ 50 ; | ln[—gr^j]! ~ 10 (6.33) 

and the exponent in fl6.32p ~ With -C 1 under the assumption of weak 

coupling and the realm of validity of the quantum master equation, we hnd that the decay 
of the power spectrum could be marginally observable if A/TT ~ 0.1 implying a suppression 
< 10% of the power spectrum for modes that re-enter the horizon near recombination, 
assuming no further corrections during the post-inhationary era. 

Therefore, although the power spectrum decays as a consequence of the interaction with 
the environmental degrees of freedom, it is likely that these corrections are of marginal 
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observational relevance, at least within the model studied here. However, this important 
observational fact notwithstanding, there is the noteworthy and fundamental aspect that the 
amplitude of the perturbation does not freeze out but decays after crossing the Hubble radius. 
These results conhrm in a non-perturbative manner previous perturbative analysis (l^. 
but also points out that not only the power spectrum does not freeze-out after “horizon 
crossing” but that the time dependence is associated with a violation of scale invariance 
even when in absence of interactions the power spectrum is exactly scale invariant. 


VII. DISCUSSION 

There are several aspects of the method and the results that merit discussion: 

• What is being re-summed?: The quantum master equation provides a non- 
perturbative resummation as is explicit in the hnal result for the power spectrum 
fl6.30l) . the question is what type of contributions are being re-summed. Tracing out 
the sub-Hubble degrees of freedom is manifest in the correlation function fl3.7ll3.8p of 
the held "0 shown in £g. O- These correlation functions dehne the one-loop self-energy 
corrections to the held y. The perturbative analysis of section fllVD shows how these 
correlation functions enter in the time evolution of the density matrix and lead to 
the production of single particles and correlated pairs and the quantum entanglement 
between the quanta of the system (x) and those of the (traced-over) environment {'ip). 
The Sudakov-type double logarithms originate from several different contributions: i) 
the logarithmic behavior of the one-loop self-energy is the same as in Minkowski space 
time because these arise from conformally coupled massless helds that act as proxies 
for degrees of freedom that remain sub-Hubble all throughout inhation. ii) the scale 
factor oc l/p in the interaction vertex, iii) the growing mode functions inside the ker¬ 
nel which lead to the factor 1 /rf in fl6.26p . which is canceled by the decaying mode 
(see eqns. (16.23116.25116.27p h thus leaving hnally the l/rj from the scale factor at the 
vertex which enhances the long time limit logarithmically (proportional to the number 
of e-folds). Thus, in summary, the quantum master equation furnishes a non- 
perturbative resummation of the secular Sudakov-type double logarithms 
of the one-loop self-energy of the y field. When the mode functions of the y 
held become “super-horizon” the quantum master equation effectively describes quan¬ 
tum entanglement between these super-horizon modes and the sub-horizon degrees of 
freedom that have been traced over, in agreement with the results of ref, fl^ . 

• Generality: We have considered the interaction oc yV'^ where y describe the degrees 
of freedom that are of cosmological relevance and become “super-horizon”, namely 
the “system” and '0 describes the sub-horizon degrees of freedom that are integrated 
out and is considered as a “bath” or environment. This interaction is, in fact, more 
general and describes various relevant cases: for example consider the case in which 
the original (unsealed) scalar held (j) features a quartic self-interaction gef)'^ during 
slow roll inhation when (j) develops a nearly constant expectation value (in slow 
roll). Writing (j){x,r]) = $ -|- 6 (j){x,r]) (here we assumed slow roll and neglected 
the time dependence of <h) and writing the spatial Fourier transform of 6 (j){x, rj) as 
6 (j)^{q,ri) + 6(j)<:{q,ri) where 6 (j)<:{q,ri) describes modes that are deep inside the Hub¬ 
ble radius all throughout inhation in terms of the mode functions (15.2p . one obtains 
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the vertex [g^) 6 (j)^{q,r]) 6 (j)^{k,r]) — q,g) which is the of form where 

X = a{g)S(j)-^, V’ = a{g) 6 (j)<. In this interpretation tracing over the (50< degrees of 
freedom to obtain the reduced density matrix is akin to a Wilsonian coarse graining 
procedure of integrating out short wavelength fluctuations leading to an effective field 
theory for long-wavelength fluctuationsjs^. A Yukawa interaction with fermionic fields 
is also described by the cubic vertex considered here. Indeed the fermionic fields are 
expanded in mode functions of the form fl2.16l) with the index i/ ~ 1/2 for a massive 
field with irtf H with technical differences associated with the spinorial contribu¬ 
tions to the loop corrections, which however yield a similar logarithmic contribution to 
the self-energy after a quadratic renormalization of the mass since the loop correction 
is similar to that in Minkowski space time. 


Minimally coupled scalars secular and infrared enhancements: If the inflaton 
(“system”) scalar field couples to other (environmental) scalar fields that are minimally 
coupled to gravity and with masses M H these “environmental” helds feature 
quantum fluctuations that are infrared and secularly enhanced when their wavelengths 
become super-Hubblej^, 0, 0 [Il|, [H, 15, 18, 2^. 

Although in this article we focused on influence of sub-Hubble correlations, the impor¬ 
tance of the infrared and secular enhancements when the environmental scalar held is 
minimally coupled to gravity merits a discussion. 

To begin with, consider the case when the environmental held is minimally coupled 
and massless in this case = 3/2 and 

p-ikrj 


u{k,r]) = 


ys ^ k/ ' 


(7.1) 


then it is clear that the correlation function K[q, g, g'] in fid.lip features logarithmic 
infrared divergences in the integration regions fc~0; |/c-|-g|~0. A non-vanishing 
mass for the held regulates the infrared, consider <C H from which it follows 
that for a minimally coupled environmental held 


A 


A = 


Ml 

3H^ 


(7.2) 


In ref. [22| it is shown that the logarithmic infrared divergences of the (self-energy) 


kernel are manifest as poles in A. The calculation of the kernel K[q,g,g'] follows the 
same steps described in detail in ref. (see the appendix in the second reference). 
The regions of integration fc~0; |/c-|-^~0 are isolated and an infrared cutoh p is 
introduced in these regions, within which the small argument expansion u{p, g) oc 
is used and the integration yields poles in A. Outside these infrared regions it is safe 
to take z/p = 3/2, the details are available in the second reference in |22|. Using the 
results of this reference we hnd to leading and next to leading order in A 

K[q, g, g'] = Kjnlq, g, g'] + AT^cb, V, v'] , (7-3) 


where 


KiR[q,g,g'] = 
Kcc[q,V,v'] = 


8 n{gg'yl‘^ A 


[qgg] 


nA 


Stt^ 
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The first line corresponds to the infrared enhancement of a minimally coupled, nearly 
massless scalar held, arising from the infrared contribution of super Hubble modes of 
the environmental scalar held in the integral in fl3.1ip . The second line is recognized as 
the kernel for conformally coupled massless helds. Inspection of the diherent contribu¬ 
tions available in ref.jl^ reveals the former arises from the regions /c~0; |fc-|-^~0 
in the integral fl3.1ip whereas the latter contribution arises precisely from the terms 
in the mode functions u{p] rf). 

With this kernel we now must obtain the non-local contribution to the quantum master 
equation from p, p'], in particular the coefficients of the respective terms in 

fl6.20l) . namely 

Jr]0 V 


Following the arguments leading to flb.lSp . the corresponding integrals are dominated 
by the region —qr]' <C 1, therefore focusing on this region and in the long time limit it 
follows that the hrst line in fl7.4p yields to leading order in A 


KiR[q,V,v'] ^ 


{q^prj 


/\2A 


47r^g^(pp')2A 


(7.6) 


With the super-Hubble behavior of g± given by fl6.25p we hnd to leading order in A 


Lo v 
r 9 -{q,v') 


'Vo 


7]' 


KiR[q,'n,'n'] dp' ~ 
KiR[q,p,p'] dp' ~ 


q3/2p2 I27i‘^(^—qpy A 
1 ln[p/p*] 

g3/2^2 12ti‘2 A 


(7.7) 


It becomes clear that the coefficients of the Q, P terms in the non-local contribution 
to the quantum master equation now feature a much stronger secular contribution in 
the long time limit p —)• 0 and are infrared enhanced by the pole in A. 


In this case we expect that the corrections to the power spectrum of the “system” 
held will feature a stronger decay as discussed in the second reference in jl^. This 
case merits further analysis and will be relegated to future study, as our goal in this 
article is to explore the impact of environmental degrees of freedom with sub-Hubble 
correlations all throughout inhation. 


In refs. it was found that self-interactions of curvature perturbations and the inter¬ 
action between a massless minimally coupled scalar held and curvature perturbations 
lead to time dependent logarithmic corrections (and powers of logarithms) to the power 
spectrum of the C, (related to the curvature) variable, conhrming previous results in 


ref.[l3|. Although we did not study the curvature perturbation and we considered 


the inhaton-like scalar coupled to a massless conformal scalar held, our conclusions 
support in a non-perturbative manner the perturbative results of ref. ji^ in that the 
power spectrum is indeed time dependent (as well as scale dependent) after “horizon 
crossing”. Furthermore, the quantum master equation provides a non-perturbative 
resummation of secular self-energy corrections (loops) and describes the asymptotic 
long time behavior of correlations well after “horizon crossing”. 
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VIII. CONCLUSIONS AND FURTHER QUESTIONS. 


We studied the dynamics of an effective field theory of two interacting scalar helds during 
inflation: a minimally coupled inflaton-like held 0 which is taken as the system and another 
scalar held ip as an environmental quantum held that is integrated out. We obtained the 
reduced density matrix for the held 0 by tracing out the degrees of freedom of the scalar held 
<p, the reduced density matrix obeys a quantum master equation which we obtained up to 
one loop in the correlation functions of the p held. The quantum master equation describes 
the decay of the vacuum state, the production of particles and correlated pairs and quan¬ 
tum entanglement between the huctuations of the 0 and p. When the huctuations of the 0 
held become super-horizon, the quantum master equation describes quantum entanglement 
between the super-Hubble degrees of freedom of the system and sub-Hubble degrees of free¬ 
dom of the environment. Renormalization aspects emerge naturally in this formulation. The 
quantum master equation provides an ehective non-perturbative description of the dynamics 
whose solution is a resummation of self-energy corrections (loops). Our main goal in this 
article is to study the ehect of environmental degrees of freedom that remain sub-Hubble 
all throughout inhation upon the power spectrum of super-Hubble huctuations of the inha- 
ton. For this purpose we considered the environmental held to be a conformally coupled 
massless held, as a proxy for degrees of freedom whose mode functions correspond to huctua¬ 
tions that remain sub-Hubble all throughout inhation and are not amplihed. In this case the 
quantum master equation provides a non-perturbative resummation of secular Sudakov-type 
double logarithms in the asymptotic long time limit. From the quantum master equation we 
obtain the time evolution of the power spectrum for super-Hubble huctuations. Even when 
the non-interacting theory features a scale invariant power spectrum, the non-perturbative 
resummation of environmental correlations (loops) leads to a breakdown of scale invariance, 
and a decay of the power spectrum at long time. Super-Hubble inhaton huctuations do 
not freeze out but decay upon “horizon crossing”. However for weak coupling between the 
inhaton and the environmental degrees of freedom, if inhation lasts for only ~ 60 e-folds, 
the corrections to the power spectrum would be just marginally relevant for cosmological 
observations. This important aspect notwithstanding, the ehective held theory description 
based on the reduced density matrix and its quantum master equation furnishes a powerful 
non-perturbative framework to study the impact of sub-Hubble degrees of freedom upon the 
quantum correlations of inhationary perturbations of cosmological relevance today. 


Further questions: The results obtained above within a simple model of interactions 
between inhaton huctuations and other degrees of freedom that remain sub-Hubble sug¬ 
gest several possible avenues of study. For example in the theory of non-Gaussianity the 


bi-spectrum is a result of a cubic self-interaction of curvature huctuations [27|, l28|. In the 


local limit for a squeezed conhguration of the momenta of the three helds, two of the mo¬ 
menta are large and one becomes small, when the small momentum (the shortest side in 
the triangle) becomes super-Hubble and the other two are sub-Hubble, the situation is akin 
to the case studied in this article. In this case the bi-spectrum describes quantum entan¬ 
glement between sub and super-Hubble degrees of freedom. Tracing over the sub-Hubble 
huctuations leads to a quadratic ehective action for the super-Hubble huctuations much in 
the same way as the ehective action for inhaton huctuations considered here. Therefore, 
it would be interesting to obtain the corresponding quantum master equation by tracing 
over the sub-Hubble curvature huctuations and study the asymptotic long time evolution 
of the power spectrum to understand if it indeed freezes out or if there is a correction such 
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as the decay as found here in the simpler scenario. We have commented on the case when 
the environmental helds are nearly massless and minimally coupled to gravity highlighting 
the enhanced infrared and secular contributions as discussed in ref. ji^. Pursuing a deeper 
understanding of this case with the quantum master equation is certainly of interest within 
this context because curvature fluctuations are described by massless and minimally coupled 
scalar helds, albeit with derivative interactions. Another relevant case to consider is that 
with the inhaton Yukawa coupled to fermionic degrees of freedom with masses <C H. As 
discussed above the mode functions for these degrees of freedom are very similar to those 
of a conformally coupled scalar held (with the ensuing spinorial structure). In the standard 
model there is a large number of fermionic degrees of freedom, which suggests a large N 
(number of fermionic species) expansion which furnishes yet another non-perturbative re¬ 
summation scheme. Tracing over these fermionic degrees of freedom yield loop self-energies 
for the inhaton huctuations, and the reduced density matrix and its concomitant quantum 
master equation can be obtained following the methods and steps highlighted here. For 
large N one can resort to the case in which the Yukawa coupling oc 1/ y/N so that the one 
loop self energy becomes exact in the large N limit. The results of some of these studies will 
be reported elsewhere. 
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